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We consider a non-coherent communication system with M 
transmit and N receive antennas in a block Rayleigh fiat fading 
channel with coherence interval of T symbol periods. We use 
the following complex baseband notation 



IB-poM. T . 1 M - |. SNH * 10 ne 



(1) 



where 5 is the T x Af transmitted matrix, X is the T x N 
received matrix, and H and W are the unknown (to both trans- 
mitter and receiver) M x N and T x N matrices of the i.i.d. 
fading coefficients and additive noise terms from CA/(0, 1). The 
transmitted symbols are also assumed to be power constrained, 

Er.,i:"- 1 E{i*. m | 3 } = p. 

The capacity of non-coherent systems was studied in [1], 
where it was shown that at high SNR or when the coherence 
interval is much greater than the number of transmit antennas, 
capacity can be achieved by using a constellation of unitary ma- 
trices. However, at low SNR. or for small values of T (e.g., 
T as L), unitary constellations are no longer optimal. The ex- 
act expression for the pairwise error probability of unitary con- 
stellations and the ChernofT upper bound are also given in [I]. 
However, these expressions appear to be intractable for the gen- 
eral case. Therefore, inspired by Stein's lemma [3], we propose 
the use of Kullback-Lcibler [3] distance between distributions to 
approximate the exponential decay rate of pairwise error proba- 
bilities. The design problem will then be to seek constellations 
that have the largest minimum KL distance between the received 
distributions assigned to their elements (maxim in code design). 
The resulting constellations coincide with the unitary designs at 
very high SNR or very low rales. But for high rate codes or at 
low SNR, different signal sets are obtained which show better 
probability of error performance. 

The KL distance between = p(X|Si) and pj — p(X\Sj), 
where S» and Sj are two different T x M matrices from the 
constellation, can be calculated as 

V(jk\\p s ) = Nu {(I T + SiSFXIr + SjSf) -1 } - NT 
- N In d e t { ( It + Si S? ) ( I r + Sj S? ) - 1 } . 

(2) 

Examining (2) shows that the KL distance consists of two parts. 
For example, for M sb 1, (2) reduces to 



(3) 



in which, the first part is due to having different magnitudes (ly- 
ing on different spheres in ), and the second part is due to 
the angle between the points (lying on different one-dimensional 
subspaces of C^). This decoupling property suggests partition- 
ing the signal space into Subsets of concentric spheres and using 
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Fig. 1 : Frame error rate comparison of 16-point real constellations 



only intrasubset and intersubsei KL distances in the maximin 
problem. 

The unitary space- Li me codes of [I] have exponential encod- 
ing and decoding complexity. In [2] suboptimaj unitary con- 
stellations with very low encoding and decoding complexity are 
proposed. In this work, we use the low complexity designs of [2] 
as the subsets of our multilevel constellation. The resulting 16- 
point real constellation forT = 2 and SNR « lOdB is a two level 
constellation, with 4 points on the inner circle and 12 points on 
the outer circle. The error rate performance of this constella- 
tion is simulated for different values of N and compared with 
the corresponding constellation proposed in [2]. The results are 
shown in Fig. I. As we see, by using multilevel constellations 
(which use contributions from both parts of the KL distance) 
instead of single level unitary constellations (in which the first 
part of the KL distance is always zero), the performance of the 
non-coherent systems can be significantly improved. Moreover, 
the new constellations have a decoding complexity similar to the 
low complexity designs of [2]. 
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